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Abstract 

The purpose of this paper is to establish the convergence in law of the sequence of "midpoint" 
Riemann sums for a stochastic process of the form f'(W), where W is a Gaussian process whose 
covariance function satisfies some technical conditions. As a consequence we derive a change-of- 
variable formula in law with a second order correction term which is an Ito integral of f"(W) with 
respect to a Gaussian martingale independent of W. The proof of the convergence in law is based on 
the techniques of Malliavin calculus and uses a central limit theorem for g-fold Skorohod integrals, 
which is a multidimensional extension of a result proved by Nourdin and Nualart in \E\. The 
results proved in this paper are generalizations of previous work by Swanson QT] and Nourdin and 
Reveillac [7], who found a similar formula for two particular types of bifractional Brownian motion. 
We provide three examples of Gaussian processes W that meet the necessary covariance bounds. 
The first one is the bifractional Brownian motion with parameters H < 1/2, HK — 1/4. The others 
are Gaussian processes recently studied by Swanson 9 , 10] in connection with the fluctuation of 
empirical quantiles of independent Brownian motion. In the first example the Gaussian martingale 
is a Brownian motion and expressions are given for the other examples. 

1 Introduction 

The aim of this paper is to obtain a change-of-variable formula in distribution for a class of Gaussian 
stochastic processes W = {Wt,t > 0} under certain conditions on the covariance function. These 
conditions are in the form of upper bounds on the covariance of process increments. For example, the 
variance on the increment on an interval of length s is bounded by C y/s, and the covariance between 
the increments in the intervals [t— s,t] and [r — s, r] is bounded by 



s 2 \t-r\- a (r- s)- + s 2 \t 



if < 2s < r < t and \t — r\> 2s, where 1 < a < § and a + /3 = §. 

For this process and a suitable function / we study the behavior of the "midpoint" Riemann sum 



The limit of this sum as n tends to infinity is the Stratonovich midpoint integral, denoted by L f (W s )°dW s . 
We show that the couple of processes {(Wj, $ n (i)) ,t > 0} converges in distribution in the Skorohod 
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space (B[0,oo)) 2 to {(W t , $(*)), t > 0}, where 

*(*) = f(W t ) - f(W ) -\J f"(W s )dB s 

and B — {B t ,t > 0} is a Gaussian martingale independent of W with variance r)(t), depending on the 
covariance properties of W. This limit theorem can be reformulated by saying that the following Ito 
formula in distribution holds 

/(W t ) 4 /(Wo) + / /'(W S )W S + i /"* /"(W.)dB,. (1) 
Jo * Jo 



The above mentioned convergence is proven by showing the stable convergence of a d— dimensional 
vector ($ n (£i), . . . , $ n (id)) and a tightness argument. To show the convergence in law of the finite di- 
mensional distributions, we show first, using the techniques of Malliavin calculus, that $ n (f) is asymp- 
totically equivalent to a sequence of iterated Skorohod integrals involving /"(Wt). We then apply our 
d— dimensional version of the central limit theorem for multiple Skorohod integrals proved by Nourdin 
and Nualart in [5]. 

Recent papers by Swanson [11], Nourdin and Reveillac [7, an d Burdzy and Swanson [2] presented 
results comparable to ([l} for a specific stochastic process. In [11], a change-of- variable form was found 
for a process equivalent to the bifractional Brownian motion with parameters H = K = 1/2, arising as 
the solution to the one-dimensional stochastic heat equation with an additive space-time white noise. 
This result was proven mostly by martingale methods. In [2] and [7], the respective authors considered 
fractional Brownian motion with Hurst parameter 1/4. In [2], the authors covered integrands of the 
form f(t,Wt), which can be applied to fBm on [e, oo). The authors of [7] proved a change-of-variable 
formula that holds on [0,oo) in the sense of marginal distributions. The proof in [7] uses Malliavin 
calculus; several similar methods were used in the present paper. More recently, [B] studied the case of 
fractional Brownian motion with H = 1/6. In that paper, weak convergence was proven in the Skorohod 
space, and the Riemann sums are based on the trapezoidal approximation. 

It happens that the conditions on the process W are satisfied by a bifractional Brownian motion 
with parameters H < 1/2, HK = 1/4. In this case rj(t) = Ct and the process B is a Brownian motion. 
This includes both cases studied in [7] and [11], and extends to a larger class of processes. For another 
example, we consider a class of centered Gaussian processes with twice-differentiable covariance function 
of the form 

E[W r W t ] = r<j) ( - ) , t>r, 



where <fi is a bounded function on [l,oo) such that 

4>'{x) = 



ip{x) 



and ip is bounded, differentiable and |^'(x)| < C(x — 1) _5 - This class of Gaussian processes includes 
the process arising as the limit of the median of a system of independent Brownian motions studied by 
Swanson in [9]. For this process, 



4>{x) — \fx arctan 



1 



y/x^T 



It is surprising to remark that in this case rj(t) — Ct 2 . This is related to the fact that the variance of 
the increments of W on the interval [t — s,t] behaves as Cy/s, when s is small, although the variance 
of W{t) behaves as Ct. Our third example is another Gaussian process studied by Swanson in [TU] . 
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This process also arises from the empirical quantiles of a system of independent Brownian motions. Let 
B = {B(t),t > 0} be a Brownian motion, where B(0) is a random variable with density / £ C°° . Given 
certain growth conditions on /, Swanson proves there is a Gaussian process F = {F(t),t > 0} with 
covariance given by 

E[F(r)F(t)] = P (r,t) = F W ^(r),B(t) < g (t)) - a' 
^ V ^ «(?(r),r) u(q(t),t) 

where a € (0, 1) and q(t) are defined by ¥(B(t) < q(t)) = a. It is shown that this family of processes 
satisfies the required conditions, where r](t) is determined by / and a. 

The outline of this paper is as follows: In Section 2, we introduce the basic environment, and 
recall some aspects of Malliavin calculus that will be used. In Section 3, a multi-dimensional version 
of a central limit theorem that appears in [5] is given. In Section 4, the theorem is applied to prove 
convergence of Section 5 discusses three examples of suitable process families. Finally, Section 6 

contains proofs of three of the longer lemmas from Section 4. Most of the notation in this paper follows 
that of [5]. 



2 Preliminaries and notation 

Let W = {W(t),t > 0} be a centered Gaussian process defined on a probability space (O, F, P) with 
continuous covariance function 

E[W(t)W(s)} = R(t,s). 

We will always assume that F is the a— algebra generated by W . Let 8 denote the set of step functions 
on [0, T] for T > 0; and let f) be the Hilbert space defined as the closure of £ with respect to the scalar 
product 

(llCi])!^^ = R(t,s). 

The mapping l[o,t] ^ W(t) can be extended to a linear isometry between fj and the Gaussian space 
spanned by W. We denote this isometry by h i— > W(h). In this way, {W(h),h £ Sj} is an isonormal 
Gaussian process. For integers q > 1, let Sj® q denote the q th tensor product of We use i} 09 to denote 
the symmetric tensor product. 

For integers q > 1, let *K q be the q th Wiener chaos of W, that is, the closed linear subspace of L 2 (Q) 
generated by the random variables {H q (W(h)),h £ ft, \\h\\fi = 1}, where H q (x) is the q th Hermite 
polynomial, defined as 

For q > 1, it is known that the map 

I q (h**) = H q (W{h)) (2) 

provides an isometry between the symmetric product space i^® 9 (equipped with the modified norm 
^7=f|| • \\sj<s>q) and T-L q . By convention, Wo = K an d Io(x) = x. 

2.1 Elements of Malliavin Calculus 

Following is a brief description of some identities that will be used in the paper. The reader may refer 
to [S] for a brief survey, or to [5] for detailed coverage of this topic. Let S be the set of all smooth and 
cylindrical random variables of the form F — g(W(4>i), . . . , W((f) n )), where n > 1; g : R™ — > K is an 
infinitely diffcrentiable function with compact support, and <pi £ Sj. The Malliavin derivative of F with 
respect to W is the element of L 2 (f2,.f)) defined as 

71 Pi 

DF = Y,Q^(W(<t>i),...,W(4> n ))4>i. 

i—1 



4 



In particular, DW(h) = h. By iteration, for any integer q > 1 we can define the q th derivative D q F, 
which is an element of L 2 (Q,Sj Gq ). For example, if F = g{W(t)), then D 2 F = g"{W{t))l® 2 . 

For any integer q > 1 and real number p > 1, let W' p denote the closure of S with respect to the 
norm || • \\jjg.p defined as 

||F||^ P =E[|FH+^E[||^||^]. 

i=l 

We denote by S the Skorohod integral, which is defined as the adjoint of the operator D. This 
operator is also referred to as the divergence operator in [8]. A random element u G L 2 (Q,Sj) belongs 
to the domain of 5, Dom S, if and only if, 



\E[(DF,u)^]\<c u ^E[F 2 } 

for any F G ED 1 ' 2 , where c u is a constant which depends only on u. If u G Dom 5, then the random 
variable S(u) G L 2 (fl) is defined for all F G D 1 ' 2 by the duality relationship, 



E [FS(u)] = E [(DF, 



u 



This is sometimes called the Malliavin integration by parts formula. We iteratively define the multiple 
Skorohod integral for q > 1 as S(S q ~ 1 (u)), with 6°(u) — u. For this definition we have, 



E[FS q (u)} =E[(D q F,u) 



where u G Dom 5 q and F G O 9 ' 2 . Moreover, if h G S) Qq , then we have S q (h) = I q {h). 
For /, g G fj® p , the following integral multiplication formula holds: 



5 p (f)6 p (g) = iTrl( P )s 2 e- 2r (f® r g h (3) 



where <g> r is the contraction operator (see, e.g., [5], Sec. 2). 

We will use the Meyer inequality for the Skorohod integral, (see, for example Prop. 1.5.7 of [8]). Let 
p fe >p(^j®' ! ) denote the corresponding Sobolev space of f)® fe -valued random variables. Then for p > 1 
and integers k > q > 1, we have, 

Wi^Wjjk-g.p < Cfc : p||u|| D fc, P ( fl ®<j) (4) 

for all u G D fe ' p (^® fc ) and some constant Ck :P . 

The following three results will be used in the proof of Theorem 4.3. The reader may refer to [5] 
and [8] for details. 

Lemma 2.1. Let q > 1 be an integer. 

1. Assume F G D q ' 2 , u is a symmetric element of Dom S q , and (D r F, <P (u)) ^ r G L 2 (fi,^' 8l9_r_: ') 
/or all < r + j < q. Then (D r F, u)q 8r G Dom S r and 

r=0 ^ ' 

2. Suppose that u is a symmetric element ofJP~*~ ' 2 (S}® 3 ). Then we have, 

D^(u)j£( k )( J \w->(D^u). 



i=0 
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3. Let u,v be symmetric functions in I$ 2q > 2 {$j® q ). Then 

t=0 w 

In particular, 

\\S"(u)\\ 2 L2{n) =E[S"(uf] =E(-) E 

i=0 ^ ' 

Proof of 1. This is proved in [5] (see Lemma 2.1). It follows by induction from the relation F5(u) — 
S(Fu) + (DF,u)f, (see [S], Prop. 1.3.3). 

Proof of 2. This follows from repeated application of the relation DS(u) = u + 5(Du), (see [8], Prop. 
1.3.2). 

Proof of 3. This follows from repeated application of the duality property, (see [5], eq. (2.12)). □ 

3 A central limit theorem for multiple Skorohod integrals 

Let X — {X(h),h e $)} be an isonormal Gaussian process associated with a real-separable Hilbert 
space Sj, defined on a probability space (fi, F, P). We assume that F is generated by X. The purpose 
of this section is to prove a multi-dimensional version of a theorem proved in [5] (see Theorem 3.1). We 
begin by defining the notion of stable convergence. 

Definition 3.1. Assume F n is a sequence of d— dimensional random variables defined on a probability 
space (f2, F, P), and F is a d— dimensional random variable defined on (f2, Q, P), where F C Q. We say 
that F n converges stably to F as n —¥ oo, if, for any continuous and bounded function / : R d — >■ R and 
bounded, R-valued, J 7 — measurable random variable Z, we have 

lim E(f(F n )Z)=E(f(F)Z). 

n— >oo 

Theorem 3.2. Let q > 1 be an integer, and suppose that F n is a sequence of random variables in R d 
of the form F n = S q (u n ) = (o~ 9 (it^), . . . , 8 q {uf l )), for a sequence of R d — valued symmetric functions u n 
in J} 2q ' 2q ($)® q ). Suppose that the sequence F n is bounded in _L 1 (f2,io) and that: 

(a) (vP n , ^>'P = i(D at F^ e ) <8 h)fje, q converges to zero in i x (f2) for all integers 1 < j,ji < d, all integers 
1 < oi, . . . , a m , r < q — 1 such that a% + • • • + a m + r = ?/ a ^ ^ £ •f)® r . 

(7>j for each 1 < i,j < d, (u\ l ,D q F^^ q converges in L 1 (J7, j^) to a random variable sy, smc/i i/iai i/ie 
matrix £ := (sij) dxd is nonnegative definite (that is, A T SA > /or a/Z nonzero A € R^J. 

TTien i 7 ^ converges stably to a random variable in R d with conditional Gaussian law AT(0, S) given X. 

Remark 3.3. Conditions (a) and (b) mean that for g > 1, some combinations of lower-order derivative 
products are negligible. For example, for q — 2, then the following scalar products will converge to zero 
in L 1 (f2,^): 

• {u\,hi <8> h 2 ) f}m for all h 1} h 2 e ft. 

• (u % n , DF£ <8> /i)^ 82 for all h £ $3 and all j (including i = j). 



(D q - 2 u,D q - z v) 



\D q 



G 



. (ui,DF> ® DF*)^, for all 1 < k,j < d. 

Only the g^-order derivative products converge to a nontrivial random variable. Usually (see Section 
6), the term {u % n , D q Ffy flig9 has the same asymptotic behavior as (u z n , u^)^ mq ■ 

Remark 3.4. It suffices to impose condition (a) for h £ So, where So is a total subset of S)® r . 
Proof of Theorem 3.2. 

As in the 1-dimensional case considered in [5], we will use the conditional characteristic function. 
Given any hi, . . . h m £ f), we want to show that the sequence 

£ n = (Fj l ,...,F*,X(h 1 ),...,X(h rn )) 

converges in distribution to a vector (F^, . . . F*,, X(hi), . . . , X(h m )) , where, for any vector A £ R d , 
Fx, satisfies 



E{e**-*-\X{hi),...,X(h m )) =exp^--A J EAj , (5) 

where A • F n = Ylj=i Aj-^n denotes the usual scalar product in R rf , and we use this notation to avoid 
confusion with the scalar product in Sj. 

Since F n is bounded in L 1 ^,^), the sequence is tight in the sense that for any e > 0, 
there is a K > such that P [F n £ [— K, K] d ^j > 1 — e, which follows from Chebyshev inequality. 
Dropping to a subsequence if necessary, we may assume that converges in distribution to a limit 
(F* ,...F*,X(hi),...X(h m )). Let Y := g(X(h 1 ),...,X(h m )), where g £ C b °°(R m ), and consider 
^n(A) = </>(A, £ n ) := E(e iA ' F "F) for A £ M d . The convergence in law of £„ implies that for each 
1 < j < d: 

lim ^ = lim iE (F^e lX F "Y) = iE (F^^Y) , (6) 

n—taa UAj n— foo 

where convergence in distribution follows from a truncation argument applied to F^. 

On the other hand, using the duality property of the Skorohod integral and the Malliavin derivative: 

^ = ,E (FiW-Y) = ,E ((<, W (e— y) )^ 
Je^ITW^")^^^ Qe(<,#V a ' f " ®D«-t) s8 j (7) 



a=0 



By condition (a), we have that (u^, D a e lX Frl ® ^ converges to zero in L 1 (il) when a < q, 

so the sum term vanishes as n — ¥ oo, and this leaves 

d 

lim zE «, YD"e^)^ q = Jim i £ E (i\ k e iX F » «, 7^)^,) 

fc=i 

d 

= -^E(A fce lA ' F » Sfej r) 



fc=i 
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because the lower-order derivatives in D q e lX Fn also vanish by condition (a). Combining this with ([5]), 
we obtain: 

d 

iE (F£e iA - F ~Y) = - J2 A fc E (e lX - F - Skl Y) . 
fe=i 

This leads to the PDE system: 
d 



dX. 



E (e lA ' F °° \X(hi), X(h m )) =-J2 ^SkjE (e iA ' F °° \X(hi), X(h m )) 



k=i 



which has unique solution ([5]). □ 

4 Central limit theorem for the Stratonovich integral 

Suppose that W = {Wt,t > 0} is a centered Gaussian process, as in Section 2, that meets conditions 
(i) through (v), below, for any T > 0, where the constants may depend on T. 

(i) For any < s < t < T, there is a constant C\ such that 

E[{Wt-Wts) 2 ] <Cisi 

(ii) For any s > and 2s <r,t<T with \t - r\ > 2s, 

|E [(W t - Wt-.)(W r - W r -.)]\ < C lS 2 \t - r\~ a {t A r - s)- p + s 2 \t - r|"i; 
for positive constants a, /3, 7, such that 1 < a < | and a + f3 = |. 
(hi) For < i < T and < s < r < T, 



\E[W t {W r+s -2W r + W r - s )]\ < 

for some positive constant Ci- 
(iv) For any < s < t < T - s 



C2S2 if r < 2s or \t — r\ < 2s 

C 2 s 2 ((r - s)-i + \t- r|-i") if r > 2s and \t - r\ > 2s 



|E[Wi(Wi + .-W , t _,)]| < 
and for each < s < r < T, 

\E[W r (W t+s - W t _ s )]\ < 
for some positive constant C3. In addition, for t > 2s, 



C3S2 if t < 2s 

C 3 s(t-s)-5 ift>2s 



C3S2 if t <2s or \t - r\ < 2s 

C 3 s{t - s)~i + C 3 s\t - r\~i if t > 2s and \t - r\ > 2s 



\E[W s (W t -W t - s )}\<C 3 s 1 i + "<(t-2s)-< 

for some 7 > 0. 



8 



(v) Consider a uniform partition of [0, oo) with increment length 1/n. Define for integers j, k > and 
n > 1: 

[3 n (j, k) —E [(Wi±i - Wjl) (w*±i - 
Next, define 

= E /3„(2i-l,2fc-l) 2 + /3„(2j-2,2fc-2) 2 ; 
!£(«)= E ^(2j-2,2fc-l) 2 +/3„(2j-l,2fc-2) 2 . 

Then for each t > 0, 

lim 77+ (t) = ?7 + (i) and lim ?7~(t) = rf (t) 

n— >oo n— »oo 

both exist, where r] + (t) , i]~ (t) are nonnegative and nondecreasing functions. 

Consider a real- valued function / e C 9 (R), such that / and all its derivatives up to order 9 have at 
most exponential growth, that is 

fW(x) < K ie ^p{K 2 \x\ a ), xeM, a<2 

for k = 0, . . . , 9, and positive constants K\, K 2 - We will refer to this as Condition (0). 

In the following, the term C represents a generic positive constant, which may change from line to 
line. The constant C may depend on T and the constants in conditions (0) and (i) - (v) listed above. 

The results of the next lemma follow from conditions (i) and (ii) . 

Lemma 4.1. Using the notation described above, for integers < a <b and integers r,n>l, we have 
the estimate, 

b 

E \PnU,k)\ r <C(b-a+l)n- r *. 

j,k=a 

Proof. Suppose first that r — 1. Let / = {(j, k) : a < j,k < b,\k — j\ > 2, j A k > 2}, and J = {(j, k) : 
a < j,k < b, (j, k) I}. Consider the decomposition 

b 

E \M,k)\= E \MM+ E \Pn{jM- 

j,k=a U,k)el U,k)e.J 

Then by condition (ii), the first sum is bounded by 

E n-*\j-k\- a < Cn-*(b-a + l), 
U,k)ei 

and the second sum, using condition (i) and Cauchy-Schwarz, is bounded by Cn~?(b — a + 1). For the 
case r > 1, condition (i) implies \/3 n (j, k)\ < C\n~^ for all j, k. It follows that we can write, 

b b 

E \Pn(j,k)\ r <C in ~^ E \M,k)\<C(b-a + l)n-^. 

j,k—a j>k=a 

□ 
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Corollary 4.2. Using the notation of Lemma 4-1, for each integer r > 1, 

(Wn(2j - 1,2k- l)f + \f3 n (2j- l,2fc-2)| r + |/3„(2j - 2,2fc- l)| r + |/3„(2j - 2,2fc-2)f) < C 
Proof. Note that 

]T {Wn(2j -1,2k- l)\ r + \[3 n (2j -1,2k- 2)\ r + \(3 n (2j -2,2k- l)f + \(3 n (2j -2,2k- 2)f) 
= £ |/?„(j,fc)f- 

j,fe=0 

□ 

Consider a uniform partition of [0,oo) with increment length 1/n. The Stratonovich midpoint 
integral of f'(W) will be defined as the limit in distribution of the sequence (see [TT]): 

L¥J 

$n(*) := X! f'(Waj-i)(Wv - W21-2). (8) 

We introduce the following notation, as used in [51: St := l[o,t]i an d 9i := lji j+i j . 

The following is the major result of this section. 

Theorem 4.3. Let f be a real function satisfying condition (0), and let W = {Wt,t > 0} be a Gaussian 
process satisfying conditions (i) through (v). Then: 

(W t , A (w t ,f(W t ) - f(W ) - i J f"(W s ) dB t 

as n — >■ 00 in the Skorohod space (B[0, oo)) 2 , where rj{t) — i] + (t) — T)~(t) for the functions defined in 
condition (v); and B — {B t ,t > 0} is scaled Brownian motion, independent of W , and with variance 

The rest of this section consists of the proof of Theorem 4.3, and is presented in a series of lemmas. 
The proofs of Lemmas 4.4, 4.5, and 4.9, which are rather technical, are deferred to Section 6. We begin 
with an expansion of f(Wt), following the methodology used in [11] . Consider the telescoping series 

f(W t ) = f(W )+ £ [f{W^)-f{WM^)\ +/(Wi)-/(Wa L¥J ), 
j'=i 

where the sum is zero by convention if I -yj = 0. Using a Taylor series expansion of order 2, we obtain 
$»(*) = f(W t ) - /(Wo) - 2 E ( Aw 1i " AW %i 

- £ i2o(W3i) + ^ fli(W^) - (/(W t ) - /(W fL¥J ) 



Y 
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where Rq, Ri represent the third-order remainder terms in the Taylor expansion, and can be expressed 
in integral form as: 

1 r vv ii 

R {W M ) = - " (Wn -u) 2 f^(u)du; and (9) 



2 



W21 



flifWai-a) = -- / " (W^ -u) 2 / (3) («)rfu. (10) 



2 . 

By condition (0) we have for any T > that 

/(ir,) /uf^j) 



lim E sup 

n-foo <t<T 



so this term vanishes uniformly on compacts in probability (ucp), and may be neglected. Therefore, it 
is sufficient to work with the term 

A„(t) := f(W t ) - f(W ) - + Rn(t), (11) 

where 

- J! /"(Wati) (^Wl - AWl-i) ; and 
i=i 

We will first decompose the term \l/ n (i), using a Skorohod integral representation. Using @ and the 
second Hermite polynomial, one can write AW 2 (h) — 2H 2 (W(h)) + 1 = 5 2 {h® 2 ) + 1 for any h e S) 
with \\h\\^ = 1. It follows that, 



*»(*) = £ /"(W%i)<5 2 (all, - 9|U 



From Lemma 2.1, we have for random variables u, F 

F5 2 {u) = 5 2 {Fu) + 26 ((DF, u)^) + (D 2 F 7 u) , 



so we can write: 



:=F n (t) + S n (t) + C n (t). 

Hence, we have A n (i) = /(W*) - /(Wo) - § (F n (t) + £ n (t) + C n {t)) + R n (t). In the next two lemmas, 
we show that the terms B n (t), C n (t), and R n {t) converge to zero in probability as n — > 00. The proofs 
of these lemmas arc deferred to Section 6. 
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Lemma 4.4. Let < r < t < T. Using the notation defined above, 

E[(R n (t)-R n (r)) 2 ] <C 





nt 




nr 




( 


Y _ 




.Y. 


) 



n 2 



for some positive constant C, which may depend on T . It follows that for any < t < T , R n (t) 
converges to zero in probability as n — > oo . 



Lemma 4.5. Let < r < t < T. Using the above notation, there exist constants Cb,Cc such that 



(B n (t) - B n {r)Y 



(C n (t) - C n (r)Y 



<c B (^ 


nt 




nr 1 \ 


y_ 




IY\) 




nt 




nr 




y 




-Y- 


) 



n 2 : and 



_ 3 

n 2 . 



It follows that for any < t < T, B n (t) and C n (t) converge to zero in probability as n —¥ oo. 

Corollary 4.6. Let Z n (t) := R n (t) - \B n {t) - \C n {t). Then given Q<t\<t<t2<T, there exists a 
positive constant C such that 

E[\Z n (t)-Z n (ti)\ \Z n (t 2 ) - Z n (t)\] <C{t 2 -t x )i. 

Proof. By lemmas (4.4) and (4.5), 



E 



{z n {t 2 ) - z n {t x )y 



< 3E 
+ 2 

< C 



{R n (t 2 ) - Rn(h)Y 

{c n (t 2 ) - cmy 



2E 



(B n (t 2 ) - BMY 





nt 2 




nti 




( 


2 _ 




_ 2 _ 


) 



Then by Cauchy-Schwarz inequality, 

E \\Z n {t) - Z n (ti)\ \ z n{t2) - z n {t)\] < (e [{Z n (t) - ZM)?] E [(Z n (t) - ZM)Y 

< c 





nt 2 




nti 




( 


2 _ 




2 _ 


) 



This estimate implies the required bound C(t 2 — ti) 2 , see, for example pQ, p. 156. 



□ 



Next, we will develop a comparable estimate for differences of the form F n (t) — F n (r). In order to 
prove this estimate, we need a technical lemma which will be used here and also in Section 6. 

Lemma 4.7. Suppose a,b are nonnegative integers such that a + b < 9. For fixed T > and interval 
[ti,t 2 ] C [0,T], let 

9a = £ &HW^)(df^ 

Then we have for 1 < p < oo 



E[\\D b g a \\P] <C 



}<8>2 

'21-2 





nt 2 




nti 




( 


_ 2 _ 




2 _ 


) 
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Proof. We may assume t\ = with < T. For each 6 we can write 



E 



fc„ l|2 



= E 



< E 



2 / (a+b) (^ 2i ^)/( a+b )(^^ i )(e^ i , £ ^ i ) (a^-a^,^-^ 



sup 

0<s<t 



sup 



£ 11— \ , £ 2m-l 



E/ Q®2 _ a<g>2 <}®2 _ q®2 \ 
\ c, 2i-l u 2l-21 u 2m-\ u 2m-2 / 
, _ \ — — —^r- — — / s>m 



Recall that condition (0) holds for / and its first 9 derivatives, so the first two terms are bounded. For 
the last term, note that by Corollary 4.2 with r = 2, 



£ 

l,m=\ 



2£— 1 u 2£-2 : u 2m-l u 2rn-2 



f)02 



= \Pn(2£ - 1, 2m - l) 2 - p n {2l - 1, 2m - 2) 2 - n (2£ - 2, 2m - l) 2 + f3 n (2£ - 2, 2m - 2) 2 

£,m=l 

nt 2 



< C 



ii 



Lemma 4.8. For < s < t < T, write 



F n (t)-F n (s)= ^(/"(^^(df^-df^ 



Then given <ti < t <t% <T, there exists a positive constant C such that 

E [\F n (t) - F n (ti)| 2 |f n (t 2 ) - K(i)| 2 ] < C(t 2 - h) 2 
Proof. First, for each n > 1, we want to show that there is a C such that, 



E 



{F n {t 2 ) - F n {h)f 



<A 


nt 2 




nt\ 










) 




_ 2 _ 




_ 2 _ 





By the Meyer inequality (TJJ there exists a constant c 2 ^ such that 



E 



(5 2 (M„)) 4 < C2 : 4||Wn|lD2,4( fl ®2), 



where in this case, 



and 



y nt 2 j 



•*ra||j}2, 4(^«2j — E||U„||^®2 + EII-DUjjII jj(g>3 + E||Z) 2 U n ||^04. 



□ 



(12) 
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From Lemma 4.7 we have E||u„||^ 2 , E||£>u„||£ 83 , E||L> 2 u„||4 84 < C ([^\ - [^l) 2 ^ 2 , and so it 
follows that, 



E 



< C 



nt 2 



rati 



From this result, given < ii < t < i2, it follows from the Holder inequality that 

E [\F n (t) - F n (h)\ 2 \F n (t2) - F n (t)\ 2 ] < (E [\F n (t) - F^h)^ (E [\F n (t 2 ) - F n (t)\ 4 ]Y 





nt 2 




rati 


) 


( 


_ 2 _ 




_ 2 _ 





As in Corollary 4.6, this implies the required bound C(t 2 — ti) 2 . 



□ 



By Corollary 4.6 and Lemma 4.8, it follows that A„(i) = f(W t ) - f(W ) - \F n (t) + Z n (t) is 
tight, since both sequential parts F n (t), Z n {t) are tight. Further, we have that Z n (t) tends to zero in 
probability, and F n (t) is in a form suitable for Theorem 3.2. In the next lemma, we show that the 
conditions of Theorem 3.2 are satisfied by F n (t) evaluated at a finite set of points. 

Lemma 4.9. Fix = t < t\ < t 2 < • • • < t d . Set F l n = F n (ti) - ^(ij-i) for i = 1, . . .d, and let 
F n = (F^, ■ ■ ■ , F^). Then under conditions (0), and (i) - (v), F n satisfies conditions (a) and (b) of 
Theorem 3.2, and so given W , F n converges stably as n — > oo to a random variable £ = (^i, . . . , £d) 
with distribution Af(0, E), where S is a diagonal dx d matrix with entries: 



T f'iWsfnids), 

J ti-1 



where rj(t) = n + (n) — n (t) is as defined in condition (v). 

Remark 4.10. As we will see later, r)(t) is continuous, nonnegative, and nondecreasing. 

It follows from the structure of S that, given W, F n converges stably to a d-dimensional vector with 
conditionally independent components of the form 

where each Q ~ Af(0, 1). Thus, we may conclude that for each i, 

k a r nw s ) dB S 
Jti-i 

for a scaled Brownian motion B = {B tl t > 0} that is independent of W t , with E [Bf] = rj{t). 

Proof of Theorem |5 To prove Theorem 4.3, it is enough to show that for any finite set of times 
= to < h < t 2 < ■ ■ ■ < td we have 

(A n (h), A n (t 2 ) - A n (h), A n (t d ) - A n (i d _!)) A (A(ti), A(t 2 ) - A(h), A(t d ) - A(t d _i)) 
as n — > oo; and that A n (t) satisfies the tightness condition 

E [|A„(t) - A n (t!)p \A n (t 2 ) - A„(t)| 7 ] < C(t 2 - h) a (13) 
for < ti < t < t 2 < oo, 7 > 0, and a > 1. 
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For A n (t) = f{W t ) - /(Wo) - \F n (t) + Z n (t), we have shown in Lemmas 4.4 and 4.5 that 

Z n {t) = R n (t) - ~ (B n (t) + C n (t)) A 

for each < t < T, and hence Z n (ti) — Z n (ti-i) — > for each U, 1 < i < d. By Lemma 4.9, the pair 
(W, F n ) converges in law to {W,F 00 ) 1 where i 7 ^ is a d— dimensional random vector with conditional 
Gaussian law and whose covariance matrix is diagonal with entries 

„2 / ill / T ir ^2„ 



4 = / f'iWsY^ds). 
Ju-1 

It follows that, conditioned on W , each component may be expressed as an independent Gaussian 
random variable, equivalent in law to 



r j"{w s )dB S 
Jti-i 



where B — {Bt,t > 0} is a scaled Brownian motion independent of W with E [Bf~\ — r](t). Finally, 
tightness follows from Lemma 4.8 and Corollary 4.6. Theorem 4.3 is proved. □ 



5 Examples 

5.1 Bifractional Brownian Motion 

The bifractional Brownian motion is a generalization of fractional Brownian motion, first introduced 
by Houdre and Villa [3]. It is defined as a centered Gaussian process B H,K = {B H ' K (t),t > 0},with 
covariance defined by, 



E[B H ,K B H , K] = l_ {t 2H + S 2H 



K 



1 

2^ 



\t - s 



2HK 



where H £ (0,1), K £ (0, 1] (Note that the case K = 1 corresponds to fractional Brownian motion with 
Hurst parameter H). The reader may refer to [1] and its references for further discussion of properties. 

In this section, we show that the results of Section 4 are valid for bifractional Brownian motion with 
parameter values H, K such that H < 1/2 and 2HK = 1/2. In particular, this includes the end point 
cases H = 1/4, K = 1 studied in [7], and H = 1/2, K = 1/2 studied in [TT] . 

Proposition 5.1. Let ^B^' K ,t>0^ denote a bifractional Brownian motion. The covariance con- 
ditions (i) - (iv) of Section 4 are satisfied for values of < H < 1/2 and < K < 1 such that 
2HK = 1/2. 



Proof. Condition (i). 



E 



= t 2HK + ^(t - s) 2HK - [t 2H + (* - s) 2H ] K -^s 



2HK 



< 



2 K 
1 



\2H 



K 



Vi-^{t 2H + (t- 
where we used the inequality a m — b m < (a — b) m for a > b > and m < 1. 



Vt—s-±(f H + (t-sY H ) K 



2 K ' 
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Condition (ii) 
E 



= ( + r2H ] K [** + ( r ~ S ) 2H ] *"[(*" S ? H + r2H ] *+[(*" s? H + (r - s) 2H } K ) 



+ _L (| t _ r + s\ 2HK - 2\t - r\ 2HK + \t-r- s\ 2HK ) . 



This can be interpreted as the sum of a position term, ^ip{t, r, s), and a distance term, jK^p(t— r, s), 
where 

<p{t, r, s) = [t 2H + r 2H } K - [t 2H + (r- s) 2H ] * - [(t - s) 2H + r 2H ] K +[(t- s) 2H + (r - s) 2H ] K ; and 

W -r,s) = \t-r + s\ 2HK - 2\t - r\ 2HK + \t-r- s\ 2HK '. 

We begin with the position term. Note that if K = 1, then ip(t, r, s) — 0, so we may assume K < 1 and 
H > j. Assume < s < r < t, and let p:=t — r. By Fundamental Theorem of Calculus, we can write 
<p(t, t—p,s) as 

2HK f [t 2H + (t-p- 2H ] (t- P - - [(* - s) 2H + (t~p- 2H ] (t-p- O 2 "" 1 ^ 

Jo 

= f f AH 2 K{1 - K) [{t - n) 2H + {t-p- 2H ] K ~ 2 (t - n) 2 "- 1 {t-p- C) 2 ^ 1 ^ dn 
Jo Jo 

< 4H 2 K(1 - K)s 2 [(t - r) 2H + (r - s) 2H ] K ~ 2 (t - rf H -\r - s) 2 "- 1 
<Cs 2 {t-r) 2HK - 2H - l {r- S ) 2H -\ 

This implies condition (ii) for the position term taking a = i + 2H > 1 and (5 = 1 — 2H. 

Next, consider the distance term ip(t — r, s). Without loss of generality, assume r < t. Again using 
an integral representation, we have 

ip(t -r,s) = \t-r + s\ 2HK - 2\t - r\ 2HK + \t - r - s\ 2HK 

= [ S 2HK[(t-r + 2HK - 1 -(t-r-0 2HK - 1 ] <% 
Jo 

= J J 2HK{2HK - 1) [t - r + rj\ 2HK ~ 2 dn d£ 
<Cs 2 (t-r- S ) 2HK - 2 <Cs 2 \t-r\-i, 



since \t — r\ > 2s implies (t — r — s) 2 < 2 2 \t — r\ 2 . 
Condition (Hi). 



E 



1 



^\[t 2H + (r + s) 2H ] K - 2[t 2H + r 2H ] K + [t 2H + (r- s) 2 ^ h 
1 



Take first the term, (p(t,r,s). If r < 2s, then 

\[t 2H + (r + s) 2H ] K - 2[t 2H + r 2H ] K + [t 2H + (r - s) 2H ] K \ < Cs 2HK = Csi 
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based on the inequality a K — b K < (a — b) K for a > b > and K < 1. Hence, we will assume r > 2s. 
If K = 1, then H = j, and we have 



I \Jr + s — 2\/r + \Jr — s\ 



1 



o 2y/r + x 



dx 



1 



1 



1 



4 Jo Jo (r - s + x + y) 

1 9 / N — 3 

< 7 s ( r - s ) 5 ; 



o 2^/r - s + x 
dy dx 



dx 



and if K < 1, 



f 2HK[t 2H + (r + xf H ] K - 1 (r + xf^dx - f 2HK[t 2H + (r - s + xf^-^r - s + xf^dx 
Jo Jo 



< 



AH K(K - l)[t 2H + (r - s + x + y) ] (r -s + x + y) 4 "' 2 dy dx 



o Jo 



+ 



[ S f S 2H(2H - l)K[t 2H + (r - s + x + yf"]^ 1 ^ -s + x + y) 2H - 2 dy dx 
Jo Jo 



<4H 2 K(l-K)s 2 {r-s) 



2HK-2 



2H{1 - 2H)Ks 2 {r - s) 2tlK - 2 < Cs 2 (r - s)"5 



This bound for (p(t, r, s) also holds in the case \t — r\ < 2s, so the bound of Cs? is valid for this case. 
Next for the second term. Note that if \t — r\ < 2s, then 



^ (|t - r + s\ 2HK - 2\t - r\ 2HK + \t-r- s\ 2HK ) 
If \t — r\> 2s, then we have 

1 



< 2(3s)^ K < Cs*. 



y/\t-r\ + s- 2y/\t -r\ + y/\t -r\-s 



.dx 



o 2yJ\t-r\ +x Jo 2y/\t-r\ - s + x 



-dx 



= ff 

Jo Jo 



< 



(\t-r\- s + x + y 



< 



dy dx 



4(|t-r|-s)5 2\t-r\z 

using the inequality \ t _l\_ s < ]t=r\ ^ or 1^ ~ r l — This bound for tp(t — r, s) holds even in the case 
r < 2s, so the bound of Cs3 when r < 2s is verified as well. 
Condition (iv). 

For the first part, we have for alH > s, 



E 



Bf K (<f - 5f_f ) 



1 [t 2H +(t+s) 2H ] K -i,[t 2H + (t-sr] K 



2 K 



This is bounded by Cs 2 if t < 2s. On the other hand, if t > 2s, 





- f 




2 K J- 



(t + x) 2iI - x dx 



<Cs{t- s) 



2HK-1 
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For < s < r < T with t > 2s and \t - r\ > 2s, 

1 



E 



K 



< 



K 1 



2H- 



K 



+ 



2 K 

±^\r-t + s\ 2HK -±\r-t-s\ 



< Cs(t-s)-i +Cs\r-t\-i. 

If t < 2s or \t — r\ < 2s, then we have an upper bound of Cs^ by condition (i) and Cauchy-Schwarz. 
For the third bound, if t > 2s, 



E 



B H.K ( B H,K B H,K^ 



< 



< 



1 



[ 8 ™+1?H] K -± [s 2H + (t-s) 2H 



K 



^(t-s) 2HK -±(t-2 S r K 



J' HK [s 2H + (t-s + x) 2H ] K (t-s + x] 
(t - 2s + x)~^dx 



2H-1 



dx 



2 K+1 



< Cs{t - 2s)~i = Cs5+T(t - 2s)" 



for 7 = \. 



□ 



Proposition 5.2. Let B H ' K be a bifractional Brownian motion with parameters H < 1/2 and HK = 
1/4. Then Condition (v) of Section 4 holds, with the functions r/ + (t) = 2C^t andr/~(t) = 2C^t, where 

C K = [ 2 +J2 (V2m + l-2V2^+V2m-l) , 

\ m=l / 

Ck= \ 2 k + i + {V2nTT2-2V2nTTT+V2^) . 



m—1 



Proof. As in Prop. 5.1, we use the decomposition, 



2- K n-^<p(j, k, 1) + 2- K n-h(j - k, 1). 



The first task is to show that 



[nt] 



lim V n"V(i,fc,l) 2 = 0. 

n— ^oo * — ' 

j,k=l 



(14) 
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Proof of ([TJ]). We consider two cases, based on the value of H. First, assume H < i. Then 



rt, k, i) = [u + i) 2H + (k + i) 2H ] K [U + i) 2H + K 

-[J 2H + (k+l) 2H ] K + [J 2H + k 2H ] K 

= f 2HK [(j + 1) 2H + (k + x) 2H ] K ~ l (k + xf"- 1 ^ 
Jo 

- f 2HK [j 2H + (k + x) 2H ] K ~ X (k + xf^dx 
Jo 



^ AH 2 K{1 - K) [(j + y) 2H + (k + x) 2H ] K ~ 2 (k + x) 2 "' 1 (j + y) 2 ^ 1 dy dx 
o Jo 

< C k 2HK - 2H - l j 2H - 1 = ck-i- 2H j 2H -\ 

With this bound, it follows that 

[nt\ [nt\ 



-1-AH 



n n ^— ' 

j,k=l j—l k—1 

n 

which tends to zero as n — > oo because H < ^. 

Next, the case H = h. Note that this implies K = h, and we have 



\<p(j,k,l)\= ^/j + k + 2 - 2^/j + k+i + y/j + k <C(j + k)-z. 
So with this bound, 

\nt\ L"*J 
n L nt J oo „ L«*J 

y TO -3<iiy r 2 

J=l m=j + l j=l 

which tends to zero as ?i — > oo because j -2 is summable. Hence, (fH)) is proved. 

From (fT4|) . it follows that to investigate the limit behavior of r)£(t), f]~(t), it is enough to consider 

- V ^(2j - 2fc, l) 2 + ^{2j - 2k, l) 2 = - V V(2j - 2k, l) 2 ; and 
n ^-^ n ^— ' 

j,k—l j,k=l 

- V V(2j - 2fc + 1, l) 2 + ^(2j - 2fc - 1, l) 2 - - V V(2j - 2fe + 1, l) 2 ; 
n n 

j.k—l j,k=l 
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since the sums of 4>(2j —2k+l, l) 2 and tp(2j — 2k — 1, l) 2 are equal by symmetry. We start with 



1 LtJ 

- ^(2. ? -2fc,l) 



n 



= i E (V|2j - 2fe + 1| - 2^2j - 2k\ + y/\2j -2k- 

= iE 4 +iEE(^- 2i +i-2^+V2j-^i)' 

j=i j=i fe=i 



4 I M I 2 L 2 J / \ 2 

= "Mr + 4^ S E (v^t+t- 2V^+ V2^r) 

j — 1 m— 1 

4 I M I o °° / \ 2 

j=l m— 1 

2 L"rJ oo 

~4^E E (^+T-2V2^ + V2^T) , 

i=i m=j 

where the last term tends to zero since 

(V2m + 1 - 2V2^ + V2m- l) < ^ (2m - 1)~ 3 < C(2j - 1)~ 2 , 



m—j m—j 

and, 

n 



]T (2, - i)- 2 



as in oo. We therefore conclude that, 

7?+(t)= lim V (/3„(2. 7 -l,2fc-l) 2 + /3„(2j-2,2fc-2) 2 ) 

n— >oo * — ' 



3,k=l 
n— too 



lim - V V(2j - 2fc, l) 2 = 2C+t, 

n— >oo n ^-^ ' 

where 



1 / °° 

C+ = ^ I 2+ ^ (\/2m + 1 - 2V2m + V2m - lj ' 

\ m=l 



For the other term, 



- V ^(2j-2fc + l,l) 2 



4 4^E( 2 -^) 2 + i^EE(v / 2j-2fc + 2-2V2j-2fc + l-V2F 
j=i j=i fc=i 



2fc 
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Hence, by a similar computation, 



r]-(t)= lim V /3 n (2j-l,2fc-2) 2 +/3„(2j-2,2fc-l) 2 = 2C7^, 

j;k=l 



where 



C 



(2 - \/2) 2 1 



A' 



+ pf H (/ 2m + 2 - 2 ^ 2m + 1 + • 



m—1 



□ 



As a concluding remark, it is easy to show that > C K , and in general we have 7/ + (i) > ry (i). 



5.2 A Gaussian process with differentiable covariance function 

Consider the following class of Gaussian processes. Let {Ft,0 < t < T} be a mean-zero Gaussian 
process with covariance defined by, 



E[F r F t ]=r<P - , t>r 



where : [1, oo) — > M is twice-diffcrcntiablc on (1, oo) and satisfies the following: 
(0.1) Moo := sup x >! \<j){x)\ < c 0! o < oo. 
(0.2) For 1 < x < oo, 



l<*'(*)l < 



(0.3) For 1 < x < oo, 



(15) 



|0"(x)| < c^-s^-l)"!. 
where c^j, j = 0, 1, 2 are nonnegative constants. 

Proposition 5.3. The process {F t ,0 <t<T} described above satisfies Conditions (i) - (iv) of Section 



Proof. Condition (i). By Conditions (0.1) and (0.2), 



E 



{F t - F t _ s ) 2 = i0(l) + (t - s)0(l) - 2(t - s)<f> 1 + 



< 2(t - s) 

< 2(t-s) 



1 + 



t - s 



J 4>{ x ) dx 



(1) 



+ s\\ 



t-s 
s 10(1)1 



C0,l 



< 2(i 

7i Vie - 1 

< Cs^Vt - S + S||0||oo 

< Cs^, 



dx + s\\(f>\\ c 
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where the constant C depends on max jVT, ||0||oo }• 

Condition (ii). For 2s < r < t — 2s we have by the Mean Value Theorem, 



I E [F t F r — F t - s F r — F t F r - s + F t ~ s F r 



t - s 



-(r-s) U 



r — s 



t - s 



r — s 



< s sup |0" 0)1 



t t-s 



r — s 



t-s\ 2 (t-s 



ts 



r(r — s) 



< 



cVT s 

(t-r)\ 



CVT s 2 \t-r\-'i. 



Condition (Hi). By symmetry we can assume r < t. Consider the following cases: First, suppose 
2s < r < t — 2s. Then we have 



\E[F t {F r+s -2F r + F r _ s )]\ = 



(r + s) 



(r + s) 



t 



r + s 
t 



r + s 



st 



< — sup \(f> {x)\ 
r _s_l 

L r+s ' r — s \ 



2r»,[ +(r-s)4 
t t 



t 



r — s 
t 
r 



r — s 



< 



< 



2sH 2 c^ 2 
r(r — s)(r + s) \ t 

CsHi 



r—s r+s 
r + s 



t — r — s 



r(t — r) 2 

There are two possibilities, depending on the value of r. If r > |, then £ < 2, and we have a bound of 



Cs 2 ( - 

r 



(t-r)< 



< 2CVT s 2 \t-r\-i. 



on the other hand, if r < ■§, then t^- < 2 and r < t — r. Then the bound is 



Cs 2 



T 



< 2CVT s 2 



(r-s) 2 + \t - r \ 2 



t — r J \ ryjt - r 

For the case \t — r\ < 2s, assume that t = r + ks for some < k < 2. Then 
\E[F t (F r+s -2F r + F r _ s )}\ 

'tV (r + s)" 



= (iA(r,p 
= (iA(r,p 

< 3(r + s) 

< 3(r + s) 



tA(r + s) 
tV(r + s) 
tA(r + s) 

(fc + l)s 



2r<j> ( - ] + (r - s)0 



r — s 



(r + s)0(l) - 2r0 ( - ) + 2r0(l) + (r - s)0 



r — s 



(r-s)^(l) 



1 + 



r — s 



-0(1) 



< 3(r + s) 



i+ 



<t>'{x) dx 



i+ 



(fc+l)3 



dx < CVTsi. 
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For the last case, note that if tAr < 2s, then we have an upper bound of Cs 2 , since E [F s Ft] < s\\ 0|| oo . 

Condition (iv). Take first the bound for E [F t (F t +s — F t - S )]- Note that if t < 2s, then an upper bound 
of Cs? is clear, so we will assume t > 2s. We have 



\E[F t F t+s - F t F t _ s }\ 



t(t> 



t + s 



(t - a)<f> 



t 



<{t-s) sup \(j)'(x) 



t ' t- 



t-s 
t + s 



t 



< c^,i — 



t It 

- + C<pfi s 



t t-s 



t + s 



t V t + s V s 

< CsVt (t-«) _ *. 

For the case r ^ t, first assume r < t — 2s. By condition (0.2) 



|E [F r F t+s - F r F t . 



r0 



f + s 



r0 



t - s 



< 2s sup \4>'{x)\ 



2s^F c^i CVT a 
~~ y/t — r — s ~ y/t — r 



lir>t + 2s, then 



|E[F r F 4+s -F r F t _ s ]| 



(i + s) 



(i-s)0 



< t 



t-s 
dx + 2s\\4>\\ 



< 



t — s + x 
2sfc^i\/F+ s 2sc^oVT 
\Jt — t y/t — s 

< Cs(r-t)-* +Cs(t-s)-*. 

For the case t < 2s or |r — t| < 2s, the bound follows from condition (i) and Cauchy-Schwarz. 
For the third part of condition (iv), we have for t > 2s, 

E [F.F t - F s F t - s ] = s0 - s0 f 1 —^- 



< s sup |0' (x) I I 

\!-JL i] V 



f t-s 
s s 



< 



< Cs'(t- 2s)"5 
= Cs^ +7 (t-2s)- 7 



where 7=5. 



□ 



Proposition 5.4. Suppose <j>{x) satisfies conditions (0.1), (0.3) and m addition (f>(x) satisfies: 



(0.4) : 0'(x) 



x/x - 1 ' 
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where itel and tp : (l,oo) — > M is a bounded differentiable function satisfying \ip'(l + x)\ < C^x 2 
/or some positive constant C^p. Then Condition (v) of Section 4 is satisfied, with rj + (t) = C^t 2 , and 
T)~(t) = Cpt 2 for positive constants C^,Cp. 

Remark 5.5. Observe that condition (0.4) implies (0.2) but not (0.3). 
Proof. We want to show 



/3„(2j-l,2fc-l) 2 -^C pA t 2 ; 
M2j -2,2k- 2) 2 — ■* C^ >2 i 2 ; and 



53 /3„(2j-l,2fc-2) 



(16) 



(17) 



(18) 



so that Cp — Cp t \ + C/3.2, and Co — 2Cp^. We will show computations for (JT6J) , with the others being 
similar. As in Prop. 5.2, 

L^J L^J L¥J j-i 

J] n (2j -1,2k- If = ]T Pn(2j ~ 1, 2j - l) 2 + 2 E ^ ~ 1, 2* - I) 2 , 



so it is enough to show 



i=i fe=i 



lim V V/3„(2j -1,2k- l) 2 = Ci^ 2 ; and 

j=l fe=l 



Urn V ^„(2j-l,2i-l) 2 = C 2 
3=1 



t 2 . 



(19) 



(20) 



Proof of (fT9|) . For 1 < fc < j — 1 we have 



2A- 



/3„(2j-l,2fc-l) = — U -7 - 



2j 



2A- 



2j - 1 
2fc 



2fc - 1 



2j 



2k- 1 



2j -1 
2fc- 1 



2j 

2fc 



0'(x) dec 
n 1 2j-i ?i 



2fc-l 



2fc— 1 



0'(x) cfa;. 



Using the change of index j = k + m and a change of variable for the two integrals, this becomes, 



Pn(2j - 1 



' 2fc - 1)= n/ + l dy -n 



2m+l 



6' 1 



2fc-l 



dy. 



(21) 



With the decomposition of (0.4), we will address (l2~Tj) in two parts. Using the first term, we have 



p2m 




f2m+l / 








/2 m -1 


/ 2/ n 


/2m V 



2k- 1 



2 k 

n 



2k (V2m - V2m - lj - V2fc - 1 (V2m + 1 - 
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We are interested in the sum 

4k 



E ~~T (V2m - V2m - l) - V2k - 1 (V2m + 1 - V2mj 



k—1 m— 1 



(22) 



We can write 



(y2m - V2m - l) - y/2k - 1 (\/2m + 1 - V2m) 



= - V2fc - 1 ( \j2m + 1 - 2V2m + V2m - 1 ) + ( V2fc - V2fc - 1 ) ( V2ro- \/2m - 1 
Observe that 



V2k - V2fc - 1 ) ( V2m - \j2m-\ 



< 



1 



(2fc- l)(2m- 1) : 



and so 



4k 



2 L^J L^J- fc 

E E 



^ Z_, ( 2fc-l)(2m-l) 



< 



4k 



2 (m 



v 



E 2 fc - 1 

fe=i 



< 



C\og{ntf 



Therefore the contribution of this term is zero, and it follows by Cauchy-Schwarz that the only significant 
term is 



4k 



2 L^J L^J- fe 

Y E (2fc - 1) (\/2m + 1 - 2y/2m + y/2m - 1 

fe=l m=l 



4k z 



I s£ I -m 

x , / . , — . \ 2 L 2k -1 

2^ f y/2m + 1 - 2V2m + \/2m - lj ^ — 

m=l fe=l 



4k 2 f\/2m + 1 - 2V2m + V2rn 

m— 1 



T \ a Oil 



which converges as n — > oo to 



oo 

K 2 t 2 Y (V2m + 1 - 2-\/2m + V2rn - l) ' 



Next, we consider the term -^ip( x )- The contribution of this term to (|2"Tj) is 



t>2m 
>2m-l 

We can bound (|2"5|) by 

,.2m 



2A: 



2k + y 



^( 1+ 2k> dy ~- 



2m+l 



2m 



2k - 1 



2fc - 1 + y V 2fc - 1 



^ 1 + dl/. (23) 



2 k 



< 



2m-l 
1 



2k + y 



2m+l 
2m 



^, V2fc-V2fc^T / 2fc 
( S) mX)l V2fc + 2m-l + V2fc + 2 TO -1 



2fc- 1 
2fc - 1 + y 

2-ro 



^ 1 + 



y 



2k- 1 



dy 



2m+l 



V> 1 + 



2k - 1 



dy 



n 
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Since is bounded, we have 



C 



< 



c 



V2fc- l v / 2fc + 2m- 1 V2fc- 1 V2m - 1 
For £fc, m using that |^i'(ac + 1)| < Cx~i, 

r 2m+l 



^ 1 + 



2fc — 1 



dy 



< 



2m- 
2 m 



/ U \ ( U + 1 \ 

/( 1+ 2ib)-^( 1+ 2*3lJ d " 
(1 + u)cfo 



2m-l 
2»i 



+ 1 
2te— 1 



< c 



< 



v 2 dv du < 



C 



V2fc- 1 



V2m + 1 - V2m 



C 



V2fc- V2m- 1 



so that 



■Bfc. m < 



2k 



C 



< 



c 



2k + 2m - 1 V2fc - lV2m - 1 ~ \/2fc - V2m - 1 
Hence, from (|2"41 and (|2"5| . we obtain 



c 



<- y 



^ (2m- l)(2fe - 1) 



< 



Clog(n) 2 



(24) 



(25) 



so the portion represented by (1231) tends to zero as n — > oo. Since this term is not significant, it follows 
by Cauchy-Schwarz that the behavior of 

L¥J i-i 

5>„(2j-l,2fe-l) 2 
i=i fe=i 

is dominated by eq. (|2"2"1) . and we have the result (TiT)|) . with 

Ci = k 2 Y (V2m + 1 - 2^2m + V2m - l) . 



Proof o/ O- For each j, 
P n (2j-l,2j-l)* = (^0(1) -2 



2j-l_ L ( 2j \ , 2j-l 
2J-1 



-<Mi) 



0(1) + (4j - 2) (0(1) - (l + ^py 
0(1) 2 , 4(2,-1)0(1) (^J 1 



2J-1 



4(2j - l) 2 



0(1) -01 + 



2J-1, 
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Since 



0-^(l + ^l) 



< 



Y by (0.3), we see that 



E 



0(1) , 4(2j- 1)0(1) 



0(1) -0 1 + 



2i-i 



which implies only the last term is significant in the limit. Again we use (0.4) to obtain: 



0'(x) dx 



2j-l 



-« ^ ^= dx J i -^(*) ^ 

2. / 1 



V2J^r V 2 i-i/' 



hence 



4(2j - l) 2 



and taking n — > 00, 



(D-0(1 + ^- 



16k 2 (2j-1) 2 /ji 
n 2 (2j-l) In 2 



.2,2 



which gives Thus (HI is proved with C^i = 4k 2 + 2k 2 Y,Z=i (V^rn + 1 - 2v / 2m + V2m - 1) 
By similar computations, 



C/3,2 = 4k 2 + 2k 2 f V2m + 1 - 2\/2m + \/2m - 1 j ; and 

m— 1 

00 2 

C/3,3 = 4k 2 + 2k 2 Y (yZm + 2 - 2V2m + 1 + VzinJ ; 



and so 



2 

C£ = C^,i + Cp i2 = 8k 2 + 4k 2 ( V2m + 1 - 2%/^ + V2m - l) , 



C /3 — 26/3,3 — 



00 

in 2 + An 2 ^ (V 2m + 2 - 2\/2m + 1 + 



Note that > , and it follows that 77(f) = rj + (t) — r\ (t) is nonnegative, and strictly positive if 
k^O. □ 

For a particular example, we consider a mean-zero Gaussian process {F t ,t > 0}, with covariance 
given by 

E\F r F t ] = Vrt, sin -1 



This process was studied by Jason Swanson in a 2007 paper [9j , and it appears in the limit of normalized 
empirical quantiles of a system of independent Brownian motions. 
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Corollary 5.6. The process {F t ,0 <t< T} with covariance described above satisfies the conditions of 
Section 4, with r/(t) = (ct -C^t 2 , where G$ , are as given in Proposition 5.4, with k 2 = 1/4. 

Proof. Assume < r < t < T . We can write, 



/■/si,. 1 ( Jj\ = Vrt tan 1 (J ) = rt f> (~ 



where 



= I ^an- (^=) , ifx>l 

If, if x = l 

Condition (0.1) is clear by continuity and L'Hopital. Conditions (0.2) and (0.3) are easily verified 
by differentiation. For (0.4) we can write, 

1 ' r^-tan"^ 1 



6' (x) = -== + —= v - tan" 



so that k = —1/2, and 

s 1 (yfx- 1 

V\ X ) = o / f ~ tan 



2 \ y/x — 1 \v^— T / 

satisfies (0.4). □ 



5.3 Empirical quantiles of independent Brownian motions 

For our last example, we consider a family of processes studied by Jason Swanson in |10) . Like [3] , this 
Gaussian family arises from the empirical quantiles of independent Brownian motions, but this case is 
more general, and does not have a covariance representation flT5j) . 

Let B = {B(t),t > 0} be a Brownian motion with random initial position. Assume B(0) has a 
density function / £ C°°(R) such that 

sup(l + \x\ m )\f n) (x)\ < oo 
for all nonnegative integers m and n. It follows that for t > 0, B has density 

u(x,t) = / f(y)p(t,x - y) dy, 



where p(t,x) — (2irt) ^e » . For fixed a £ (0, 1), define the a-quantile q(t) by 

u{x, t) dx = a, 

i 

where we assume f(q(0)) > 0. It is proved in [10] (Theorem 1.4) that there exists a continuous, centered 
Gaussian process {F(t),t > 0} with covariance 



'(B(r) <q(r),B(t) < q(t)) 



a 



2 



E [F r F t ] = p(r, t) = " > \ J \ K ■ (27) 

u[q[r),r) u{q{t),t) 

In [10] , the properties of p are studied in detail, and we follow the notation and proof methods given 
in Section 3 of that paper. Swanson defines the following factors: 

p(r, t) = P (B(r) < q(r), B(t) < q(t)) - a 2 ; and 9(t) = (u(q(t), t))' 1 ; 
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so that p(r, t) = 6(r)6(t)p(r, t). For fixed T > and < r < t < T, the hrst partial derivatives of p are 
calculated in [in](see eqs. (3.4), (3.7)): 



Q Mr) 

— p{r. t) = q(t) J p(t-r,x- q(t)) u(x, r) dy dx 

I /•«(*) 
~ 2^(* ~ r ' < ?( r ) ~ Q(t))u(q(r),r) + u(q(r),r)q'(r) / p(t ~ r, q{r) - y) dy 



2 



— oo J — oo 



d 2 

p(t-r,x-y)—u(x,r) dx dy; (28) 



d 1 



(29) 



Lemma 5.7. Let < T , and < r < t < T . Then there exist constants Cj, i = 1, 2, 3, 4, suc/i i/iai: 
(a) 



p{r,t) 



0>) 



(c) 



(d) 



d 2 



d_ 

dt 

dt; 2 



P(r,t) 



P(r,t) 



<C x \t-r\ 



<Co\t- 



<C 3 \t-r\-2 



< d\t-r\ 



Proof. Results (a) and (c) are proved in Theorem 3.f of [10]. Bounds for (b) and (d) follow by differ- 
entiating the expressions for d r p{r,i) and d t p{r,t) given in the proof of that theorem. □ 

Proposition 5.8. Let T > 0, < s < T A 1, and s < r < t < T . Then p(r,t) satisfies conditions (i) - 
(iv) of Section 4- 

Proof. Conditions (i) and (ii) are proved in [Hi (Corollaries 3.2, 3.5 and Remark 3.6). For condition 
(hi), there are several cases to consider. 
Case 1: s < r < t — 2s. Using Lemma 5.7(a), 



|E [F t (F r+S - 2F r + F r ^ s )} \<\p(r + s, t) - p(r, t)\ + \p(r, t) - p(r - s, t)\ 



< 



p(r + x, t) 



dx 



d 



p(r + y,t) 



dy 



< 2 / Cx\t-r-x\-2 dx < Cs'- 
Jo 



Case 2: If \t — r\ < 2s, the computation is similar to Case 1, where we use the fact that 

dx = 2s5. 
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Case 3: For r, t > 2s and \t — r\ > 2s, the results follow from Lemma 5.7 (b) and (d) for r < t and 
r > t, respectively. 

Now to condition (iv). For the first part, we first assume t > 2s. Then using the above decomposi- 
tion, 

E [Ft(F t+s ~ Ft- a)] = P(t, t + a)- p(t, t - s) 

= 6(t) [9(t + s)p(t, t + s)- 9{t - s)p(t, t - s)] 
= 9{t) [p(t, t + s)A9 + 9(t - s)Ap] , 

where A9 = 9(t) - 9(t - s) and Ap = p(t, t + s)- p(t, t- s). First, note that 



\u'(q(t),t)\ 



-^u(q(t),t)q'(t) + ^u(q(t),t) 



where we used that q'(t) is bounded (see Lemma 1.1 of |10j). Since u(q(t),t) is continuous and strictly 
positive on [0, T], it follows that 9{t) is bounded and 



hence, 



For Ap we have 



\m\ = ^ff < c, 

u 2 {q(t),t) 



\A9\ < J \9'(t + x)\ dx < Cs. 



(30) 



\Ap\ = |P (B(t) < q(t), B(t + s)< q(t + a)) - P (B(t) < q(t), B(t - s) < q(t - s))\ 



q(t) r q(t+s) 



< 



< 



< 



— oo J — oo 
q(t-s) r q(t) 



q(t-s) r q(t) 



— oo J — oo 



oo •/ — oo 
q(t-s) 



p(s, x — y)u(x, t) dy dx 

p(s, x — y)u(x, t) — p(s, x — y)u(x, t — s) dy dx 
\u(x, t — s) — u(x, t)\ dx + Cs 



p(s, x — y)u(x, t — s) dy dx 
Cs 



r f 9 . yj , _ i f°° r d 2 , NJ 

/ / —u(x,r)dr dx + Cs = - / — ?u[x, r)dr 

J-oo Jt-s or 2 J^^ J t _ s dx 2 



dx + Cs 



1 />00 />00 j°t 

-2 

^ J —oo J —oo Jt—S 

When t < 2s, we write 



|/ (y)\p(r, x — y) dr dy dx + Cs < Cs. 



\E [F t (F t+s - Ft-,)] | < \ P (t, t + s)- p(t, t)\ + \p(t, t) - p(t - a ,t)\ 

) f" d 

I —p(t,t + x) dx + J —p{t + y,t) 

< Csi, 



dy 



using Lemma 5.7 and the fact that 



2 dx = 2s 2 . 



30 



For the second part of condition (iv), we consider 

\E[F r (F t+a - F t - a )]\ and |E [F s (F t - F t - a )]\ ■ 

When r < t — s (including r = s), an upper bound of Cs\t — r\~^ is proved in [10] (see Corollary 3.4 
and Remark 3.6). When r > t + 2s, or \t — r\ < 2s, the bounds follow from Lemma 5.7. □ 



The rest of this section is dedicated to verifying condition (v). We start with two useful estimates. 
As in Proposition 5.8, suppose 0<s<r<t<T. It follows from Lemma 1.1 of 10 that for some 
positive constant C, 

\q(t)-q(r)\<C(t-r). (31) 
Using this estimate and the fact that e~° — e~ b < b — a for < a < b, we obtain 



(q(t)-q(r-s))* 



e 2(t-r) _ g 2(t-r+s) 



< Cs < 1. 



(32) 



Recalling the definitions in condition (v), we can write for t € [0, T] 

ti(t)-Vn(t)= E M*-M-l) 2 + 2 ^ /3„(2fc-l,2j-l) 2 + 2 £ /3„(2fc - 2, 2j - 2) 2 
«=l k<j-l k<j—l 

-2 £ ,3„(2fc-2,2j-l) 2 -2 ^ /3„(2fc-l,2j-2) 2 . 

k<j-l k<j-l 



For the first sum, since Fe_ — Ft-i is Gaussian, we have 



By Theorem 3.7 of [10], 

Lntj 



( E [( 



F £ — F&-i 



3 



F i — F t—\ 



£ (f£ - — > - / (u(q(s), s))' 2 ds 

a i o 



in L 2 as n — > oo. For the second sum, assume 1 < fc < j, and we study the term 



,'2k 2j\ (2k - 1 2j 

p n 2fc-l,2j-l )=p — ,^ -p ,^ 

n n I \ n n 



2k 2j - 1 



2 / 

2j-l 



9'(r)p ( r, ^ + 0(r)9 r p (r, - 
n ) \ n 



-P 
dr 



2k - 1 2j - 1 



n n 



9\r)p r. 



2J-1 



Q(r)d r p r 



We can write this as 



2j 



, , r) ( ^rP ( r, — ) - d r p ( r, — — - 

n J 1 2k-i \ \ n J \ n 



2j-l 



9{r) \d r p\r 



dr 

2i-i 



(//• 



9'(r) 



2j\ _f 2j 
— P r, — 
n / \ n 



2j-l\ _( 2j-l 
P r, 



dr. 



dr. 



(33) 
(34) 
(35) 
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The first task is to show that components (34) and (35) have a negligible contribution to r)(t). For (34), 
it follows from (|3U| that 



2j-l 



< Cn-\ 



(36) 



and using (|29|) . we have 



(r)9 r p r, 



2i-i 



cir = 



2j-l (2j-l 
r, q 



— q(r) I dr < Cn 



Hence, the contribution of (34) to the sum of /3„(2fc - 1, 2j - l) 2 is bounded by C [n 2 ^ ■ n 2 < Cn 1 . 
We can write component (35) as 



?'(r) 



2j 



2.7 



2J-1 



2i-i 



p r 



2J-1 



dr. 



Using (|29]>. we have for each r € f 2 ^-, ^] 



P r, — - p [r, 



<Cn-*(2j -2k-l)~*. 



Then, using (j3"6) and ([3D]), we have (35) bounded by 



C 



71 ~ 2 (2j - 2k - 1)"2 +71" 



Hence, the contribution of (35) to the sum of /3„(2fc — 1, 2j — l) 2 is bounded by 

Cn" 2 [ n ~ 1{ ?j - 2k - ^ + n ~ 2 ] ^ Cn ' 1 - 

j=l k=l 

We now turn to component (33). By (|29|) . 



d ( 2j\ 1 (2j (2j , 



or V n I 2 \ n 



To simplify notation, define 



ipn(j,r) = e 

By ([3"T]). we have for the interval I 2 k = [? k ~ 1 2fel 



(g(t)-"( r » 

2(i-r) 



n 7 n 



) ■n ' 



sup 



((g(¥)-g(0) ) < gX2j-2fc + l) 



2(2i 

V n 



This implies that inf (2j, r) , r 6 I 2 k} > 
is close to unity. We can write, 



0(r) ( d,./5 I r, — ) - d r p I r 



2J-1 



<ir 



hence, when j, k are small compared to n, \t/)\ 

dr (37) 



1 



2V27T JilLzl 2j 



2j-l 
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! ; (l-^ n (2j-l,r)) [ * ' \dr (3S) 



2^ J 21^1 " ' \Jli-r ./lEI_ r 



1 ^ 2 JZ M2j,r) ~ M2j ~ 1,r) "~ 



For component (38), by the above estimate for inf {ip n (2j, r),r e 7 2 fc} we have 

sup |1 - tp(2j,r)\ < Cn- 1 (2j - 2k + 1) < 1, 

hence (38) is bounded by 

Cn'i {2j - 2k + 1) (y2.7 - 2fc + 1 - 2^2.?' - 2k + yj2j - 2fc - l) . 
Given e > 0, we can find an M > f such that 



2 

^ (V2to + f - 2\/2m + V2m - f ) < e. 



m=M 

\2 



The contribution of (38) to the sum of (3 n (2k — 1, 2j — l) 2 is thus bounded by 

{2™y x V 9 2 - V sup (1 - V„(2j, r)) 2 f v/2j - 2fc + 1 - 2^2j - 2A: + v^j - 2fc - l) 

l_TrJ j-M-1 2 

< Cn" 1 -2k + 1- 2yj2j - 2k + y/2j - 2k - l) 

i=i fe=i 

L^J i-i , 

+ Cn- 1 Y H Cn" 1 (2j - 2fc + 1) ( a/ 2 ! -2k + 1- 2^2j - 2k + ^2j - 2k - 1J 

j'=l k=j-M 

i=i j=i 

which is less than Ce as n — > oo, since is bounded. 

For (39), by we have sup {\ip n (2j, r) — ip n (2j — 1, r) | ,r € hk} < Cn _1 , hence (39) is bounded by 
Cn~2 [2j — 2k — 1)~ 5 . Therefore the contribution of the term including (39) to the sum of (3 n (2k — 
1, 2j — l) 2 is bounded by 

Cn" 3 ^ J^(2j - 2ft - I)" 1 < Cn~ 2 log(nt), 
i=i fe=i 

because #(t) is bounded. 

It follows that the sum of (3 n (2k — 1, 2j — l) 2 is dominated by (33), and the significant term in (33) 
is given by (37). Hence, it is enough to consider 

2 raflj 



tvk z — ' \ n 

j<k-i 



Y 02 {^) ( V 72 .? - 2fc + 1 - 2 v/2j - 2fc + v^j - 2fc - l) ' 
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Using the change of index j = k + m, this is 



rnr 



E 



2j 



^ (y2m + 1 - 2v/2m + V2rn - l) . 



m— 1 



Taking n — > oo, this behaves like 



where 



By similar computation, 



- f 6 2 {s) ds, 

T JO 

oo 

a = ^ (^/2m + 1 - 2\/2m + \/2m - l) ' 



rn—1 



V /3„(2fc - 2, 2j - 2) 2 — ► - f e\s) ds, 



k<j-i 



b /"* 

V (3 n (2k-2,2j- If — > — / 2 (s) ds, and 

V /3„(2fc-l,2j-2) 2 ^^ / 2 ( S )ci S , 



fe<j-l 



where, 



00 , ,2 

bi=Y^ W 2m + 2 - 2V2m + 1 + V2^J , 

m— 1 

b 2 = (V2m - 2y/2m - 1 + \/2m - 2) . 



We have proved the following result: 

Proposition 5.9. Under the above assumptions, p(r,t) satisfies condition (v) of Section 4, where 

2 + 4a - 2b! - 2b 2 /"* , . , 

V{t) = / {u(q(s),s)) ds. 

k Jo 

The coefficient 2 + 4a — 2bi — 2b 2 is approximately 1.3437, while u(q(t),t) depends on / and a. 

6 Proof of the technical Lemmas 

We begin with two technical lemmas. The first is a version of Corollary 4.2 with disjoint intervals. 
Lemma 6.1. For < t < h < t 2 < t 3 < T, 



hm £ £ (df^-dfl^df^-dfl,) 

n— >oo ^— ' ^— ' \ — ^ — — — — ^ — — n — / $ 

j= L^aJ+ifc=L^aj+i 



= 0. 



Proof. We may assume to = and t\ = t 2 - Observe that 

= f3 n (2j -1,2k- l) 2 - p n (2j -1,2k- 2) 2 - n (2j -2,2k- l) 2 + /3„(2j - 2, 2fc 
Therefore, it is enough to show that, 

[nt 2 \ L™*3j 

E E /3™Ufc) 2 <Cn- £ 
3=0 fc=L» it 2j+i 

for some e > 0. We can decompose the sum in (j40]) as: 

\nt 3 \ \nt 3 \ Vnt 2 \-1 \nt 3 \ 

E /3n(0,fc) 2 + E /3n(N 2 J,fc) 2 + E E Pnti,k) 2 . 

k=[nt 2 ]+l k=[nt 2 ]+l 3=1 fc=|n* 2 J+l 

By condition (iv), for some 7 > we have 



[nt 3 J [nt 3 \ 

E ^(0,fc) 2 < sup |/3„(0,fc)| E 
fc=Lnt 2 J+i i<J<L«t 3 J fc=L«t 2 J+i 

[_nt 3 \ 



fc=|ni 2 J+2 

By condition (ii), for some 1 < a < |, 

E /3 n ([nt 2 \,k) 2 < P n ([nt 2 \, [nt 2 \ + l) 2 + Cn" 1 E M[nt 2 \,k) 

k=[nt 2 \+l k=[nt 2 \+2 

L"* 3 J 

^Cn^+CVi" 1 E (k - lnt 2 \)~ a < Cn- 1 , 

k=[nt 2 \+l 

and again by condition (ii), for /3 = | — a, 



[nt 2 \-l lnt a \ 

E E 

j=l fe=Lnt 2 J+l 



\nt 2 \-l \nt 3 \ 



E E Pn&kf^cn- 1 e E [(*-Ln*2j)- a r* +(*-ir s 



j=l fe=[nt 2 J+l 



[nt 3 \ 



[nt 2 \ 



[nt 2 \ 



< Cn" 1 I E * _0 J IE i^J +Cn" 1 E (N»J -i) 



hence the sum is bounded by Cn 6 for e = min {/3,7, 5}. 
Lemma 6.2. For < t <T and integer j > 1, 



< CrT* 



for a positive constant C which depends on T . 
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Proof. By conditions (i) and (ii), we have for j > 1 and t > 0, 

\nt\-l 



£ t — £|n*J 



oo 

< C^rT* (\j - k\- a A 1) + 0(n~2) < Cn^ 



(41) 
□ 



6.1 Proof of Lemma 4.4 

By the Lagrange theorem for the Taylor expansion remainder, the terms Rn(W2j ), -Ri (\V 2j-2 ) can be 
expressed in integral form: 

1 r w 2i 

Ro(W?i) = - " (W2_ L -u) 2 f^ ) (u)du; and 

2 JW 2j -i 



After a change of variables, we obtain 



1 [W*±L 



' (W2 3 -2 -u) 2 f {3 \u)du. 



and 



Define 



and 



1 r 1 

Rn(W2 3 ) = -{Wji -W2 3 -i) 3 / v 2 (vWu-x +(1- vWv ) dv\ 

2 « « J 

1 r 1 

fll (W^jj-a ) = -(W2J-2 - W2 3 -l) 3 / ^ 2 /' (3) (^2 j -l + (1 - «)P72J-2 ) <fo. 
n 2 » » Jo 

1 f 1 

G (2j) = / ^ 2 / (3) (^^ + (l-v)W2_ 1 )dv; 



1 f 1 

i(2j-2) = - / ^f^faW^-i +(1-101^2,-2)^. 



We may assume r = 0. Define AW± — Wt+i — W± . We want to show that 

2-1 



E 



< C 



y 



(42) 



This part of the proof was inspired by a computation in [7] (see Lemma 4.2). Consider the Hermite 
polynomial identity x 3 = H 3 (x) + SH^x). We use the map S q (h® q ) = q\H q {W{h)) (see © in Sec. 
2), for h £ Sj with ||/i||i5 = 1. For each j, let Wj := \\AW±\\^, and note that condition (i) implies 

Wj < Cn~^ for all j. Then 



A^f /AlfA / AW, \ id 
- H* I * 1 + WA i > - * 3 



'Si 



3 - 



5 -V + 36 
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so that 



AW? = i<5 3 (0f 3 ) + ^<5(5 X ). 

re Z 7T J » 



It follows that we can write, 

Gg(2j)AWl-i - G 1 (2j - 2)£Wl-, 

= G (2j)<5 3 (9|^) - G 1 (2j - 2)S 3 (df^) 

+ 3wl j G (2j)S(d2_ 1 _i) - 3u^-_ 1 Gi(2j - 2)^(9^2 ). 

It is enough to verify the individual inequalities 



E 



E 



E 



53 G (2j)S 3 (d^ 1 ) 

3=1 

53 G!(2j-2)5 3 (^ 

3 = 1 
3 = 1 



< c 



T 



< c 



< c 



nt 
~2 



nt 

T 



_ 3 

n 2 . 



and 



E 



53 wl j _ 1 G 1 {2j-2)5{d M ^ 

3 = 1 



< c 



nt 
~2 



We will show (|43J and (j45j) . with (|44j| and (gBj) essentially similar 
Proof of (|43|) . Using ([3]) and the duality property, 

E 



LV 



53 Go(2j> 3 (a|^) 

3=1 



(43) 



(44) 



(45) 



(46) 



= E E 

3,fc=l 
L^J 3 

j,fe=l r=0 



\r=0 



Go(2j)G (2fe) ( Y,S 6 - 2r (df^T ® SfjT) 

E (^-^(Go^Go^fc)),^ 3 ^®^ 



For integers r > 0, we have 



D r G (2j) = D r J \v 2 f {5) (yWv-i + (1 - t))^) 
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2 J « 2 / (3+r) (wWati + (1 - (wfti + (! - w ) £ i r ) ( 47 ) 



By product rule and (|47|) we have 



E 



D e - 2r (G (2j)G (2fc)) , dflr ® df, 3 -7 



(48) 



a+f)=6-2r 



sup 

0<D,U)<1 



1 /•! 



X 



f (a) (^2j-i + (1 - vW*i-*)& ) (wW2 k -i + (1 - w)^) 
+ (1 - t/)e|°) ® {we^ + (1 - W )4?) , 0fjr ® df^ 



du dt/;. 



Notice that by condition (0), E sup |/ (3+r) (£)| P 

< oo, where the supremum is taken over the random 

variables {£ = vW Sl + (1 — v)W S2 , < t> < 1, < s±, S2 < T} . From Lemma 6.2, for integers a and b 
with a + 6 = 6 — 2r, we have the following estimate 



( (uef^ + (1 - v)efj) ® (ttfefL + (1 - ™)4? ) , 0f?T ® df 3 ^ 
It follows that if r ^ 0, then by Lemma 4.1, Equation (|48|l. and Equation (|4U)) 

> 6 - r (G (2j)G (2fc)) , dfgT ® df 3 ^ 



dvdw < Gn~ (3 " r) . 

(49) 



E 



j,k=l 



4=>0S- 



< G 



«-3 



which satisfies (|42|) because § — 3 < — |. On the other hand, if r = 0, then 



£ Cn' 3 < C 

j,k=l 



nt 
~2 



and we are done with (1431) . 

Proof of (|4"5j) . Proceeding along the same lines as above, 



E 



2- 




). 


= E 



E2 2 



G (2j)G (2fc) {<5 2 (Wx (8>9 2 .:-i) + (d^-i.d^k-x) } 



< Gn" 1 £ E 

j,k=l 

+ Gn" 1 £ E 

3,k=l 



E sup IGoWl' 



d 2j-l , 8 2k-! 



£ E\(D a Gn(2i)D b G n (2k),d 2 (dy-i ® 3 



.a+&=2 
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By Lemma 4.1 we have an estimate for the second term: 



2j-l , U 2k-1 



< c 


nt 









_ 3 

n 2 . 



Then the first term has the same estimate as (I48| when r = 2, which proves (|4"5)l and the lemma 



6.2 Proof of Lemma 4.5 

As in Lemma 4.4, we may assume r = 0. Start with B n (t). Define 



7, 



.(*) ■■= J2f i3) ( w ^{ £ ^ 



7 ^2^-1 u 2j-2 



J'=l 



LttJ 

d ij-i — 8 23-2. 

j=1 



so that B n {t) = 2<5( 7 „(i)). By Lemma 2.1, we have ||tf(7*(i))ll£»(n) < E ll7n(i)||| + E||£»7 n (t)||» ea . We 
can write 

LttJ 



Il7n(*)||«= E / (3) (W / ^)/ (3) (W A -^)(^,a2 ti -0^ 



X ( £ 2fe — 1 , p 2k-l — p 2k-2 ) ( 8 23-I — 8 23-2 , 8 2k-! — p 2k-2 
\ n n re / £ \ n n n i» / 



< sup 

0<s<t 



/^(W s ) SUP ( £ 2j-l , G> 2j-1 — ^ 2j-2 ) (t? 2j-l — C? 2j-2 , C> 2fc-1 — & 2k-2 



< c 




1 




n 2 




y_ 





Note that E [sup 0<;s<t |/( 3 )(W S )| 2 ] = C by condition (0), and by Lemma 6.2, \(et,dy-i - d v-i 
C2"nT^ for all j, t. By Corollary 4.2 we know, 

^ y ( p 2j-l — 8 23-2 , 8 2k-! — d lk-2 
.!■>■' I 

Hence, it follows that E||7„(*)||| < C [f\ ra^n"* < C [f\ n~i. Next, 

Dj n (t) = ^2 f^( W*j-i ) l^£ 2j~i , 8 2,-1 — 8 23-2 \ (e vj-i ® (d zj-i - 8 23-- 



< 
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and this implies 

\\D ln (t)\\%^ < sup f^(W s )\ 



0<s<t 



^ ' /g 2j-l , d 2j-l — 8 23-2 ^ 2fc-l , 8 2k-! — d'2 



3,k=l 



£ 2j-l <g) ^G> 2j-1 — ^ 2j-2 , £ 2fc-l ® ^ 2fc-l — G> 2fc-2 



< sup 

0<s<t 



f (A \W t ) 



sup ( £ 23^1 , 8 23-1 — 821 



0<s,r<t 



3,k=l 



x sup |(e s ,£r) fl | ^ ( ^ 2,-1 — 8 23-2 ,8 2k-! —8 



By condition (0), E [sup 0<s<t |/^ 4 ^ (Ws)|] is bounded, and sup 0<s r<t | (e r ,£ s )^ | is bounded. Hence, it 
can be seen that E||Z?7„(i)||~^ 82 gives the same estimate as 7«(t). 

For C n (t), using condition (0) and the identity a 2 — b 2 = (a — b)(a + b), we can write 



E[C n (t) 2 ] <E 
By Lemma 6.2, 



sup \f^(W s )f 

0<s<t 



£ 2j-l , 8 23-1 — 8 





( 




sup 




\l<j<!f 



■2j 2 



£ 2 j - 1 , 8 23-1 — 8 23 -2 \ S ' ( £ 23-1 , 8 23-1 + 8 23 - 
^ 3 = 1 



< C^n 2 t and by condition (iv), 



\ 



E 

i=i 



£ 2j-l , l r 2j-2 2j 



j=2 



nt 
T 



Hence it follows that E [C„(i) 2 ] < C |_^J n 2 for some constant C, and the lemma is proved. 



6.3 Proof of Lemma 4.9 

For i = 1, . . . , d, set 

i=L^J+i 

and recall that = <5 2 (u^). As in Remark 3.3, we want to show: 
Condition (a). For each 1 < i < d, the following converge to zero in 

(a.l) (u l n , hi ® h 2 ) ^ m for all hi, /12 S .fj of the form £ T (see Remark 3.4). 

(a.2) (<, £>i™ ® ft) ^ for each l<j<d andheSj. 

(a.3) <<, CF„ fc ) fl82 for each 1 < j, fc < d. 
Condition (b). 

(b.l) «,£ 2 F^ 2 -^OinL 1 ifi^j. 

(b.2) (u % n , D 2 Fn)^® 2 converges in L 1 to a random variable of the form 



FL = c t r(W s ) 2 V (ds). 
Jti-1 
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The proofs of (a.l) and (a. 2) are essentially the same as given in [3] (see Theorem 5.2) but the proof 
of (a. 3) is new. 

Proof of (a.l). We may assume i = 1. Let hi ® h 2 — e s ® s T £ f)® 2 for some values s,r G [0, i]. Then 

(E> fa)^ = X! /"O^Hizi) ( 9 2^i _ 92j^i,e s ^ (d^j-i - 822-2,6^ ; 
3=1 



so that 



(u] l ,hi®h 2 ) s . (S2 < sup \f"(W s )\ sup sup ((92,-1 - 823-2, e s ) (^2,-1 - 823-2, e T 

' 0<s<t i<j<[2£l.j 0<s<ti \ " B 'i5l ~| l\ » 

It follows from condition (iii) that for fixed r > 



^ ( 8 21-1 — 8 23-2 , S T 



3=1 



3 = 1 



3=2 



< Cn"5 



(50) 



and Lemma 6.2 implies, 



so that 



sup sup 

l<j<L2*ij 0<s<t 



G? 2j-1 — C? 2j-2 , £ ; 



E 



0/12)^^ I) <Ct in x ^0. 



Proof of (a. 2). As in (a.l), assume i = 1. Using the same technique as in (a.l), we can write DF^ ® h 
as DPI ® e T for some r G [0,T]. By Lemma 2.1, DF* = D5 2 (ui) = 8 2 {Dul) + 5(vP n ), which gives 



{u l n ,DFl®e T )^ 2 = (u 1 nl 8 2 {Dui)®e T )^ 2 + (u 1 n ,5{u^)®e T ) 



i5»2 



For the first term, 



E 



(<,<S 2 (L><)®£ T ) fl82 U V E f"(W2i-i ) (821-1 -8 2,-2. 5 2 (DuD) -82,-2.6. 



< 2E 


sup 


E 


sup 




0<s<ti 




_1<^<I_2±LJ 



8t,5 2 (Dui) 



^ ' ( 8 21-1 — 8 21-2 , £ T 



By (J50J), the sum has estimate Cn 2 ; and for the second term we can take 



8 L ,5\Dui)) < sup \\d A U \\S 2 (Dui)U. 

It follows from condition (i) that — Cn~*. This leaves the ||5 2 (D3t^)|j^ term. By the Meyer 

inequality for a process taking values in S), 



E [\\S 2 (Dui)\\%] < Cl E\\Dui\\%» 3 +c 2 E||DVj| 2 04 + c 3 E||D 3 <|| 2 « 5 , 



(51) 
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so that by Lemma 4.7, E [||<5 2 (-D«)||^] < C, and we have 



E 



(ul5 2 (Dui)®s T ). 



< Cn—i. 



Then similarly, 



< 2 



sup |/"(W S )| sup (dt_,8(u 3 n )\ (d ze-i — 8 21-2 , st 



0<s<ti 



Similar to the above case, for each 1 < t < I 



9t,S{u 3 n ) 



< E 



|9a|UII«5(<)|U 



hence with (l5Cfl) we have 



E 



< Cii^i (E||<|| 6 ®2 +E||D<|| i5 i83) < Cn~i 



Proof of (a. 3). For this term we consider the product (u l n ,DFi <£> DF^) Lemma 6.1 shows that 
scalar products of this kind are small in absolute value when the time intervals are disjoint, therefore 
it is enough to consider the worst case (m* , DF^ ® DF^) ~® 2 > and assume ti = i. We have 



E 



£=1 L 



£=1 

<c£e[ 



! n 



8 21-1 — 8 21-2 , DF n 



l,2£-2 2l ] ,DF n 



Using the decomposition DF 7 \ = 8 2 (Du i n ) + <5(m*), the above summand expands into four terms: 



(1) 


( 821-1 


— 821^ 




>. 








(2) 


( 821-1 


— 82c- 


-_2 ) S 2 (Du 1 n ) 


>. 






I 




(3) 


( 822-1 


— 821- 


-Ji,fi( u l)) a 




(' 








(4) 


( 821-1 


— 821^ 


z2.^(K)) 
1 




(' 
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We will show computations for the terms (1) and (4) only, with the others similar. For (1) we have 



^ E E [ - d2_i_2, 8 2 {Du 1 Sj /l [2< -2 M ] ,S 2 (Du 1 n ) 



C E { d2t - 1 -d2 t -2,6 2 (f< 3 )(VK2 m - 1 ) £ 2 m -i (dfLi 

,S 2 (,f (3) (^ : - 1 ) £2m: - 1 - 

s 2 (/ (3) (^)(^-a|^)) 



0) 



< C sup E 

i<fe<L¥j v 

X E / C> 2<;-1 — (9 21-2 , 

£,m,m' = 1 



£ 2m -1 



Lj2£-2 21 j , £ 2rn'-l 



By Lemmas 2.1 and 4.7, the Skorohod integral term is bounded by Cn 2, and we use conditions (hi) 
and (iv) for the scalar products to obtain an estimate of the form 

Cn' 2 E ((2m -I)"* +|2*-2m|-*) ((2* - 2)~* + \2£ - 2ro'|-*) < Cn"i . 

£ 7 m,m f — 1 

For term (4), we have by a computation similar to the proof of Lemma 4.7, 



E 



5 (f (3) fVF 2fc -i) (dgfc-i - 92*-2^ 



and by conditions (i) and (ii) we have 

Cn~i ^ ^ 2^-1 — C92£~2 , C?2m-1 — (3 2 m- 

<Cn-i E (|2^-2m|- Q ) (|2£-2ra'|- a 

£,m : m' — 1 



21] , <9 2 to'-1 — C?2 



Proof o/ f&.ij. By Lemma 2.1, we can expand D 2 F n as follows: 

«, £> 2 i^ 82 = «, «5 2 0D 2 <))^ 2 + 4 «, W))^, + 2 «, 



(52) 



Without loss of generality, we may assume that u l n is defined on [0,£i] and F£ is defined on [t\,t2\ 
for ti < t2, so that the sums are over 



< = E /"(^) - s&) ; ^d < = E /"(ww) (all, - dfljj ■ 



E 

m=L^J+l 
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First term 



E 



«,<5 2 0D 2 <))^ 2 



■ E 



< E 



< E 



. ! nt 2 j 

E f"(w^) (dfL - dlh) ,s 2 I E /^(w^je^ ® (afL, 

Vm=L^j+l 



sup \f"(W s )\ 

0<s<t 



sup 

0<s<t 



E 



EE 

. £ m 



-®2 

' 2m-l ) ^2^-1 U 2£-2 



S 2 (f (4) (Wi m -i)(dT*_i-dft 



I I I rtt 2 

2 J L 2 



SUp||(5 2 (34)||L2(f2) E E 



2 I \ 2 



f=l m=l 

First we need an estimate for the 5 2 (gA) term, where in the notation of Lemma 4.7, 

34 :=/ (4) (W^) (df^-df^). 

By Lemma 2.1, \\6 2 (g 4 )\\ L 2^ < c 1 E||3 4 || fl 02 + c 2 E||£>54|| n ®s + c 3 E\\D 2 g 4 \\^ m , and so ||<5 2 (.g 4 )||L2(f2) < 
Cn~5 for each [^J < m < [^J- We can write, 



E 



«,<5 2 0D 2 <))^ 2 <Cn"i £ 



l,m=l 



£ 2m-l , (? 2l-l ^ — ^£ 2m-l , 9 21- 



We need an estimate for the double sum. We have by condition (iii), 



E 



=i l 



£ 2,ii-l , 8 -21-1 \ — (^£ 2 m -l , G? 2l-2 J 



< sup 



£ 2m— 1 , 1 j- 2£ — 2 2l j y l ^ ^ l (^E 2m-l , p 2£-l — $2£-2 



<Cn _ * E (K-m|"5 + {I- \)~\\ A 1 

e,m=i 

I " f 2 I 



£=1 n=l 



This provides an upper bound for the double sum, hence the first term of (|52p is 0(n 2 ). Note that 
in the above estimate the double sum is taken over 1 < £, m < I • It follows that this estimate also 

holds for the case i = j, that is, E (<,<5 2 (L> 2 <)) 



i5®2 



Second Term 

Using t\ < t2 as above, 



E 



= E 



E /"(^)( 



2j-l U 2j-2 



,5 



E / (3) (^)^®(ari-^ 
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E 



< CE 



f'(Wv-i ) Un-i , - 823-2 \ (dv-i -82,-2,82,-1 -0a-a^ fl^/WpPa-i) (d a *-i -02*- a ) 



j fc 



sup |/"(W S )| 

0<s<t 



sup 



sup||J(ff3)|U»(n)) 



j=0 fc=0 



where in this case, gs corresponds to the term including f^ 3 '(Wt)- It follows from Lemma 6.2 that 
sup I (e s ,8k/ n )^ I < Cn~2- and the double sum is bounded by Cn? by Corollary 4.2. This leaves an 
estimate for H^^IIl^si)- By Lemma 2.1, PCs^ll^n) < Ci||g 3 |U + c 2 1 1 -0^3 II fl® 2 ■ For this case, 



Ml < 



hence ||<?3||j5 < Cn 4. Similarly, 



sup |/ (3) (W fl )| S 

,Q<s<t 



p2k — l — p2k-2 



< Cn -3 , 



sup |/ (4) (W 5 )| 

0<s<t 



sup ||e s ||)5 


82k-! — 82,-2 


0<s<t 





hence the second term is 0(n 4 ). As in the first term, the double sum estimate shows that this result 
also holds for (<, 8(DF l n ))^ 2 . 

Third Term 
We can write 

«,<)^ 2 < sup \f{w s )\ 2 £ £ (a?, 2 , -dfl 2 ,df^ -8fL 2 ) 



and it follows from Lemma 6.1 that 



SjX2 



< Cn £ , for some e > 0. 



Proof of (b.2). As in case (b.l), this has the expansion ([52]) . From remarks in the proof of (b.l), the 
first two terms have the same estimate as the i ^ j case, hence only the term (u l n , u^)^ 82 is significant. 

Third Term 

Assume for the summation terms that the indices run over L%ij + 1 < j, k < L^J ■ We have 

j,k " » 



) ^ 2fc~l L/ 2fc-2 



Expanding the product, observe that, 

- agjli - ) = /3„(2j - 1, 2fc - l) 2 - /3„(2j - 1, 2fc - 2) 2 

- /3„(2j - 2,2*; - l) 2 + /3 n (2j - 2, 2fc - 2) 2 , 

where wi) is as defined for condition (v). For each n, define discrete measures on {1, 2, . . . }® 2 by 

00 

M + := ^ £„(2j - 1, 2k - l) 2 + /3„(2j - 2,2fc - 2) 2 ^ fe ; 



M - := 2 /3 n (2j - 1, 2k - 2) 2 + /3„(2j -2,2k- 1) 2 5, 



fe- 
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where in this case Sjk denotes the Kronecker delta. In the following, we show only n+, with n n being 
similar. It follows from condition (v) that for each t > 0, 



/i+([0,t] 2 ) := lim u n 





nt 




nt 




( 


_y_ 




_~2_ 


) 



= lim V f3 n (2j - 1,2k -l) 2 + (3 n (2j- 2,2k -2) 2 =n+(t). 

n ' 

Moreover, if < s < t then 

) + £ £ /3„(2j-l,2fc-l) 2 +/3 n (2i-2,2fc-2) 2 

which converges to /i + ([0, s] 2 ) because the disjoint sum vanishes by Lemma 6.1. Hence, we can conclude 
that n n converges weakly to the measure given by u + ([0, s] x [0, t}) = n+(s A t). It follows by continuity 
of f"(W t ) and Portmanteau Theorem that 

/" (/3„(2j - 1, 2k - If + p n (2j -2,2k- 2) 2 ) 





ns 1 






ns 












) = Mn ( 


_T_ 







= / / // (W a )/ // (W r tl )la< t l u <t/i+(d*,d«) 

JR 2 



converges to 



/ f"(W s fn + (ds). 
Jo 

Combining this result with a similar integral defined for u~ , we have for t > 0, 
lim V/"(Wti)f(H't 1 )te 2 , 

n— >00 z ' n n \ — ~ — 

j,k=l 

ft 



2j-2_ J ^2fc-l " 2fc-2 



= f f{W s ) u+{ds)- f f"{W s ) u-{ds)= f f"(W 8 ) V (ds) 
Jo Jo Jo 



where we define rj(t) = n+(t) — n (t). It follows that on the subinterval {ti-i,U\ we have the result 
in as n — > 00. □ 



!"(W s ) 2 n{d S 

Jti-X 
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